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Abstract 

We develop a semiclassical theory for the long range excitation of plasmon 
resonances in atomic clusters, based on the doorway hypothesis. The effect 
of the width of the plasmon resonance is fully taken into account. As an 
application we study plasmon excitation in small Sodium clusters, in collisions 
with electrons and protons. 

I. INTRODUCTION 



Metal atomic clusters have been studied intensively for many years to explore the differ- 
ence between atomic-like and bulk material behavior of the nano-meter metal particles [|TJ . 
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Intermediate size clusters, being between molecules and bulk materials, have shown very 
complex energy spectra and relaxation dynamics. This complexity results from chaotic be- 
havior of multi electron energy levels and their interaction with multiple vibrational modes 
of the cluster atoms. Experimental data about photo-absorption and photo-ionization spec- 
tra of mesoscopic clusters have exhibited pronounced structure related to the collective 
excitation of cluster electrons [0,3. To explain the manifestation of plasmon-like modes in 
small size atomic clusters, "ab initio" calculations have been done in the framework of the 
Functional Density Method (FDM) and the Random Phase Approximation (RPA) 0- [|J. 
Later, a better agreement between theory and experiment has been reached using Hartree- 
Fock and RPA evaluations for detailed treatment of electronic collective modes |5|||. At the 
same time, semi-empirical models have been developed for simple explanations of the ex- 
perimentally observed collective phenomena on the base of the surface and volume plasmon 
effects 

An alternative way to investigate the dynamics of excited cluster states is through clus- 
ter excitation by charged particles: protons, ions, fast electrons or charged clusters. The 
involvement of a large number of internal degrees of freedom influences strongly the cluster 
response to the external field. For example, electronic excitation may relax so rapidly among 
numerous vibrational modes, that it can change the character of a collective wave function 
during an excitation process. We will show, that the dynamics of cluster excitation can 
be described as the excitation of some envelope state of cluster eigenfunctions, even under 
conditions when the cluster spectra display a strong chaotic behavior. The scale alteration 
of the time and intensity of the cluster excitation by variation of velocity v, impact param- 
eter b and charge Z of projectile particles may be used for investigations of the dynamics of 
the cluster excitation and its energy relaxation. Analysis of the energy for different (v, Z) 
in conjunction with a description of the structure of the chaotic part of clusters' spectra 
can provide the basis for a theoretical description of excitation and decay of the collective 
excitations in mesoscopic clusters. 

In the present article we have developed a semi-classical method for the description of the 
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long-range cluster excitation in collisions with charged particles. Our theoretical approach 
takes fully into account the finite lifetime of the excited modes and their reconstruction 
into chaotic states of excited clusters. The dynamical description presented in this article 
may be considered as concrete realization of the doorway phenomena |§|9|]. Excitation of 
the doorway mode is well known for inelastic collisions of nuclear particles as well as for 
electronic transport through mesoscopic quantum dots [jTOfl . 

The paper is organized as follows: The basic semiclassical reaction theory is presented in 
Section II. In Section III we give a full account of the doorway model of plasmon resonance 
excitations. This model is then incorporated, in Section IV, into the reaction theory of 
Section II, and application of the theory to several cases are then presented in Section V. 
Finally, in Section VI, concluding remarks are made. 

II. PLASMON OSCILLATIONS AS SEMICLASSICAL COULOMB EXCITATION 

OF ATOMIC CLUSTERS 

We consider in this section the excitation of collective plasmon resonances in metal 
clusters. The probe could be any charged particle (electron, positron, proton, nucleus, 
cluster etc.). The theoretical framework of our development is a time-dependent semiclassical 
collision theory [|TTJ. The collective modes excited in the atomic clusters are treated as 
damped vibrations specified by an energy and a decay width. The novel feature of our 
theory is the treatment of the life-time of the resonances and reconstruction of the chaotic 
part of the cluster in a dynamical way. 

We show in figure 1 a schematic representation of our scattering system. The interaction 
between the projectile, with charge Ze, and the cluster's free electrons is given by 

N Ze 2 Ze 2 Ze 2 N - 

where N is the number of free electrons in the cluster, R(t) is the unit vector parallel to 
the projectile's position vector, R(t), and r, gives the position of the i-th free electron inside 
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the cluster. Above, the dipole approximation has been introduced. The projectile motion 
is described classically and, therefore, once the over-all mean field potential describing the 
projectile-cluster interaction is given, R(i) is considered to be known. The interaction, 
Eq. ([!]), can be rewritten as 

V(r,,R) = -D-E(t), (2) 
where E(t) is the classical time-dependent electric field vector, 

m = ^ am (3) 

and D = e J2i=i r i is the quantum mechanical dipole operator. 

Having specified the degrees of freedom of the system, we turn now to the description of 
the collision process. The time-dependent Schrodinger equation that describes the evolution 
of the atomic cluster is 

[H + V\ V(r 1 ,T 2 ,..,T N ,t) = ih ^ , (4) 

where H Q is the intrinsic Hamiltonian of the cluster whose spectrum is given by 

H \<Pn) = £n \<Pn) ■ (5) 

Making the expansion 

m=J2 a n(t)e- iE " t/n \<Pn), (6) 
n 

we derive the following familiar set of coupled equations (we are taking e — 0) 

ft M = ^ e -^ 4WanW (7) 

dt n 

ih = e^ n V n0 (t) a (t) + £ e^-^M Vnm{t) flm(t)| . (8) 

These equations are solved with the usual initial conditions 

a n (t -> -oo) = <5 n0 , (9) 



and the excitation probabilities are 

Pn= |a n (t^+00)| 2 . (10) 

The expansion in Eqs. (|7|||) is overwhelmingly complicated in the case of mesoscopic 
clusters, especially in the region of dynamic excitation. A convenient way to describe the 
average behavior of the dynamics is through the introduction of doorway states, a well 
known concept in the theory of nuclear reactions [§]. In the following, we give an account 
of this approach. 

III. DOORWAY STATE DESCRIPTION OF DYNAMICAL COLLECTIVE 
EXCITATION OF ATOMIC CLUSTERS 

It is well known that quantal many-body systems exhibit both collective and statistical, 
chaotic, behaviors. The interplay between these modes is clearly exhibited in the response of 
these systems to external probes. One sees usually in the spectrum a concentration of dipole 
oscillation strength at a given excitation energy which is spread over a rather large interval. 
The spread defines the width of the so-called collective excitation which is connected to the 
coupling to the more complex, chaotic degrees of freedom. 

Both in the theory of nuclear reactions and in the many body theory of atoms Q- 
|| , it has been common to give explicit reference to the collective state, although it is not 
an eigenstate of the isolated quantal system. The term 'doorway state' is used to label 
the collective state usually called Giant Resonance in the nuclear physics context. The only 
eigenstates of the systems are the complicated configurations mentioned in the introduction. 
However, due to the usually very large number of these states, it has been proven useful the 
use of the doorway state to describe the average excitation of the system. In the following, we 
develop a theoretical framework that allows the treatment of the excitation of the doorway 
state (collective state) in the collision process. Implicit reference to the statistical degrees 
of freedom is allowed through explicit consideration of the damping width of the doorway. 



5 



The starting point of our consideration is Eq. (|5|), which describes the spectrum of the 
atomic cluster. In a collision process, the cluster is originally in the ground state, \<po), and 
transitions to excited states take place through the external coupling, V. We shall assume 
that these excited states are modulated by a doorway, which could be a plasmon resonance. 
Further, to reach these 'fine structure' states, we assume that the systems has to go through 
the doorway(s) \d) HH- We formally accomplish this by expressing \ip n ) as 

\<p n ) = a n \d)+J2Pn H , (11) 
m 

where the states \m) form an orthonormal set which spans the intrinsic subspace orthogonal 
to \d). Notice that neither \d) nor |m) are eigenstates of the system. The doorway nature 
of \d) is accomplished by requiring that (ip \V\m) = 0, so that the ground state is coupled 
to the excited states through 

(ip \V\ <p n ) = V 0n = a n (ip \V\ d) = a n V od (12) 

Above, we assume that a n = (d\ip n ) depends exclusively on the energy of the state \<p n ), 
namely, ot n = a(e n ). 

It is convenient to represent the intrinsic Hamiltonian in the space spanned by the door- 
way and the subspace orthogonal to it, in the following way: 
H = \ip )e (<p Q \ +^|m)A m (m| + \d) e d (d\ + ^ { \m) A m (d\ + \d) A* m (m|} . 

m m 

(13) 

The last term in the above equation represents the spreading of \d). We emphasize that, 
although (m\d) = 0, neither \d) nor \m) are eigenstates of the system. We have also ignored 
coupling among different |m) in H . 

It is now a simple matter to derive an expression for a(e n ) and A m . Inserting Eq. (|13D 
into Eq. (§), and using that < yv|y?n >= 8n'n, we obtain 
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HSn)\= (14) 

-I _|_ v-i 



and 



T -n — £d + 



I A, 



A, 



(15) 



It is convenient for our purposes to obtain an expression for |a(e n )| which contains 
explicit reference to the difference e n — e^- For this purpose one must eliminate A m from 



Eqs. (|I4D and (|15|). This is analytically possible only for a restricted set of cases. If one uses 
the model of ref. namely a constant A m = A, and a uniform spectrum [A m ] in the sense 
A m = ms, where s is the constant spacing between adjacent levels and m = 0, ±1, ±2, 
exact expressions for e n and |a(£: n )| 2 can be derived: 

ttA 2 



cot 



7T£, r 



(16) 



and 



a\e r 



(ttA/s) CSC 2 (7T£ n /s) i + (7rA/g) 2 + {£d ^ 2 n) 



(17) 



Defining the damping width of the doorway state T d = 27rA 2 /s = 2nA 2 p, where p = 1/s is 
the constant density of the \m) states, we may finally write 



lot [e r 



2ixp 



1 d 



9 V 12 

{Sn ~ £d) + ~f + A 



(18) 



The term A 2 in the denominator is usually absorbed in the definition of the doorway energy: 



e' d = Ed + A 2 . This then results in the celebrated Breit-Wigner formula for \a(e r 
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(e n -e d - Ae d ) 2 + - 



(19) 



We may now normalize a (e n ), 



E 



a [E r 



J I a (e n ) | 2 p(s n ) de n = 1 



(20) 



The normalization above may be easily performed if we make the assumption that the 
density of the eigenstates of the system, p(e n ), equals that of the states orthogonal to the 
doorway, p = 1/s. 



Since the doorway decays into the open channels, one should add to 1^ in Eq. (@) the 
escape width, T^. This can be formally done using reaction theory techniques ||. Therefore, 



is the total width of the doorway state. 

The doorway approach described above depends on the incident energy. The faster the 
colliding probe is, the less details will be exhibited by the collision spectrum. If the energy 
is lowered, one expects that more details of the cluster structure will be revealed. Namely, 
the doorway state described above will fragment, at lower energies, owing to the coupling 
to some particular \m) states. 

The method proposed here should be particularly useful in the collisional determination 
of the complicated spectrum of large atomic clusters. This is because a detailed description 
of such large many-body system is, in general, prohibitively complicated. 

IV. PLASMON OSCILLATIONS AS DOORWAY CHANNEL FOR COULOMB 
EXCITATION OF METALLIC CLUSTERS 

Let us now consider the specific case of an electric dipole resonance. If (fo has angular 
momentum j — and the coupling V is not too strong, only excited states with j = 1 (in h 
units) are significantly populated. In this case, such states can be labeled by their energies 
(in principle a continuous label) and the z-componente of j, fi = 0, ±1; namely: n = e,fi. 
Thus, we get the coupled channels equations 




where 




(22) 




(23) 



S 



da e^) __ 1 iet h v f f \ 

at in 



(24) 



with the initial condition a n (t — > — oo) = 8 n fi. Above, we have used the fact that the dipole 
coupling does not mix the excited states. This is so because the dipole operator is a spherical 
tensor of rank 1 and all these states have j = 1. In this case the excitation occurs through the 
three degenerate doorways \d^ =0j ±i), with energy e d and angular momentum j = 1. These 
doorways differ only in the z-component of the angular momentum, /i = 0, ±1. We assume 
that the expansion coefficients a(e) are the same for the three doorways, independently of 
fi. The coupling matrix elements are then written 



V^, = a{e) (d M | V{t) \<p ) = a{e) V^t) . 



(25) 



Integrating Eq. (|24D over time, substituting the result in Eq. (|23|) and using Eq. (|2~5"D, we 
obtain 



da (t) 1 



dt 



-2 E Ids \a{e)\ 2 [ dt' e-"W* [V^t')]* a (t f ) . 

Jl ,, J J~ oo 



(26) 



Assuming the Breit-Wigner shape for \a(e) | 2 , the integration over e can be carried out in 
the complex plane. We find, 

dao (t) 



dt 



nE^) f dt ' e-^-^/Dit-O/n ao{l/) 
n ,. J —oo 



(27) 



This equation can be solved by introducing the auxiliary amplitudes v4 M (t), such that 
ao (t) = i+J2u Afj, (t). Eq. (|27D then reduces to the set of three coupled differential equations 



K (*) + 



v,(t) h\ d 2 d 



K (*) + 



v,(t)\ 2 
h 2 



(28) 



with initial conditions (t — ► — oo) = 0. Eq. can be handled by standard numerical 
procedures. Once theasymptotic amplitude ao (t — > oo) is calculated, the excitation proba- 
bilities are easily obtained from Eqs. C |l0]) and fl24j) (integrated over time) 



= |a eAt (t -> oo)| 2 = — 

Z7T 



dfe fet / R [V^(t)]* a (t) 



(29) 



Integrating over the impact parameter (implicit in the dependence of V(t) on the trajectory), 
we obtain the excitation cross sections 

^ = 2vrW bdbP £ ,,(b) (30) 

In Eq. (|30l), cr^ stands for the cross section for population of states with quantum number /i, 
and b m i n represents the value of b below which trajectories suffer strong inelastic transitions. 
We take b min to be an effective radius of the cluster core ~ I A/" 1 / 3 where / is the mean 
interatomic distance in the cluster and M is the number of atoms. 



V. APPLICATIONS 

Now we apply the formalism of the previous sections to the plasmon resonance in Sodium 
clusters, identified here as the doorway state. The resonance is excited in collisions with 
electron and proton projectiles. We consider the excitation in both neutral and ionized 
clusters. In the former case, one can use straight line trajectories for R(£) while in the latter 
one should use Coulomb trajectories. 

First we write the explicit expression for the coupling, V^. Using the dipole approxima- 
tion and writing Eq. (Q) in terms of the spherical components of the dipole field E^(t) and 
of the dipole operator D^, we get 

V = -J2 E,(t) ■ D,, (31) 

with 

E,(t) = E(t) ^ Y^(R(t)) , (32) 
where E(t) is the electric field strength 

m = ^y (33) 

If we use the reference frame of figure 1, where the z-axis if orthogonal to the collision 
plane, the x-axis is along the apex-line towards the projectile and the initial velocity has 
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positive y-component, the matrix element V^=o vanishes identically ||11|| . Using Wigner 



Eckart's theorem and representing the other two components with the shorthand notation 
Vfjt=±\ = V±, we obtain 

V + = \^o\\D\\d)\ m e _ im ^ (34) 
V- = Vft . (35) 

The phase 0(t) is given by 

where Z c is the charge of the projectile and G is the scattering angle. 

The reduced matrix-element (ipo 1 1 D 1 1 d) can be estimated from the energy weighted 
sum rule, applied to collective oscillations in the small sodium clusters H- H. It leads to 



(tpo 1 1 D 1 1 d) = e 



y ^ , (37) 



where K is the fraction of the total oscillator strength exhausted by the plasmon mode, 
M is the number of atoms in the clusters, and m is the electron mass. As an illustration, 
we consider a small Sodium cluster with Af = 8 atoms. In this case, we have the typical 
values: Ed ~ 2.5 eV = 0.093 a.u., = 0.1 Ed, K ~ 0.7, and the interatomic distance 
/ ~ 3 A. Since we are dealing with a peripheral process, we use as the lower limit for the 
impact parameter integration of Eq. ( |30| ) the radius of the cluster. In the present case, 

bmin Rcluster ~ 7.5 a.U.. 



A. Excitation of a neutral atomic metal cluster 

Let us consider the excitation of the plasmon resonance in a metal cluster with net charge 
equal to zero. In this case it is reasonable to use a straight line trajectory for the projectile 
motion. In the reference frame of figure 1, the trajectory reduces to a straight line parallel 
to the y-axis, so that 
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R(t) = tx + vty. (38) 
The electric field strength (Eq. Q3HD) and the phase (Eq. fl3"SD) are then given by 

m - < 39 > 

and 

<j>(t) = arctan f . (40) 

We have solved numerically Eq. (^) taking as projectiles electrons with incident energy 
E = 100 eV and protons with E = 100 keV. In figures 2a and 2b, we show the time evolution 
of the probabilities for the cluster to remain in the ground state, Po, and to be excited to 
a state with each of the possible values of /i, P\ and P_i. The latter where obtained by 
summation over the excitation energy. The results where obtained for the impact parameter 
b = 7.5 a.u. It is interesting to notice that the population probability for jj, — — 1 is much 
larger than that for \i — 1. 

In figures 3a and 3b, we give the calculated excitation spectra for the same projectiles 
and energies. We notice that the results are rather similar, despite of the different sign of 
the projectile charge in each case. The reason is that the cross section depends basically on 
the projectile velocity, which is of the same order of magnitude in the two cases. We have 
found that effect of the width T d (taken here as T d = 0.1 TiUd [0,(|,[7[]) on the cross sections is 
very small. For all practical purposes one may set 1^ = in the coupled channel calculation, 
using the doorways as eigenstates of H , and then use Eq. (|29|) for the probabilities P £fJi . 

Consistently with the probabilities of figures 2a and 2b, we find that the contribution of 
the state with \i — — 1 completely dominates the cross section. This can be understood if 
one neglects the width Y & and use first order perturbation theory to calculate the excitation 
amplitudes. In this case the amplitudes a±(t) = a^ = ±i are obtained by integrating Eq. (||) 
with ao(t) = 1. Using the explicit form of the matrix elements (Eqs. (|33|) and (|34"D), we get 



i < ipo 1 1 D 1 1 d > 



a±(t) = -- ^ " / dt' E(f) exp 
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h 



(41) 



Using Eq. ( |3~9|) and ( [TOD we find that the stationary phase points of the integrand are given 
by the condition 

hv 



£d = T 



b 2 

(42) 



b 2 + v 2 t 2 ± _ 

where t± stands for the stationary point in the integrand giving the amplitude a±. Since 
the energy Ed is positive, a+ has no stationary point. For a_, the stationary point is 



hb b 2 

t- = J - 2 . 43 

V ve d v 2 

Within the stationary phase approximation, the amplitude a+ vanishes, consistently with 
the low probability and cross section for /i = +1 in figures 2a, 2b, 3a and 3b. 

It is important to point out that the valitity of the stationary phase approximation 
requires that the phase in Eq. (f|l]) varies rapidly within the collision time, r c . The opposite 
situation would be the sudden limit where the doorway energy is negligible, i.e. Ed « h/r c . 
In this case, it can be easily seen that a_ ~ a + and a ~ a,. 



B. Excitation of a charged metal cluster 

We now consider the excitation of the plasmon mode in a ionized metal cluster with net 
charge Z c e, in collisions with the same projectiles of the previous section. All the parameters 
of the metal cluster are the same as in the previous subsection. In this case, R(t) is a 
Rutherford trajectory. We first consider the case where the Coulomb field is repulsive (Z 
and Z c have the same sign). We introduce the Sommerfeld parameter 

ZZ c e 2 

rj = , (44) 

the half-distance of closest approach in a head-on collision, 

ZZ c e 2 
a = ^2E~ 

and the excentricity, 



(45) 
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(46) 



It is conveniente to introduce a new variable w, defined by its relation with the time, 



t — — [e sinh w + w] 



(47) 



^From the equations of motion under the action of the Coulomb field, it is easy to show that 
the trajectory on the :n/-plane, specified by the coordinates R and (see figure 1), is given 
in terms of w as 



R(w) = a [e coshw + 1] 

/ a \/l — e 2 sinh w 



[w = arcsm 



R(w) 



(48) 



Using the variable w in Eq. (^) we get 
1 



4»+ ,, , 



y y (t(w)) i t _% 
v,(t(w)) h\ £d 2 d 



WW) + 



hf(w) 



o, 

(49) 



where primes stand for derivation with respect to w and 

r , , dw v 
f( w ) = -JT 



dt R(w) 
^From the amplitudes (w), we obtain 



(50) 



a (w) = l + J2 A » H 



(51) 



and 



a efi (w) = ^- [a{e)}* 



jet( w ')/n [ Vlt (v/)]* a Q (w') 



dw' 

1W) 



For attractive Coulomb interactions, it is necessary to make the changes: 



u r -i 

t — - [e sinh w — w\ 



(52) 



R(w) = a [e cosh w — 1] 



(53) 
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the other equations remaining the same as for the repulsive Coulomb field. 

We have performed calculations for the excitation of the plasmon mode in the same 
cluster and for the same projectiles as in the previous sub-section. However, here we consider 
an ionized Sodium cluster. In this case, the projectile's trajectory is no longer a straight 
line but a hyperbola. To assess the importance of this effect we have calculated the total 
plasmon excitation cross sections (integrated over excitation energy and summed over ji) as 
a function of the collision energy and of the cluster's net charge. In this simple example, we 
have disregarded the effects that the cluster charge might have on its structure. 

The results are shown in figures 4a, for electrons, and in 4b, for protons. In collisions 
with electrons, we notice that the cross section decreases with the collision energy in the 
energy range considered. Although the deviations from results using a straight line trajectory 
(curve with Z c = 0) are very small for high energies, Coulomb effects on the trajectory play 
an increasingly important role as the collision energy decreases. The situation for proton 
projectiles is somewhat different for the collision energies considered in figure 4b. First, we 
notice that the cross section reaches a maximum at E ps 25 keV. The second difference is 
that the cross section is not affected by the cluster's net charge. In this case the curves for 
0, 1, ...,8 units of charge cannot be distinguished in the figure. This result is not surprising. 
If we estimate the maximal deflection angle O involved in the graph: for maximal charge 
( Z c = 8), minimal impact parameter (b = b min = 7.5 a.u.) and minimal energy (E = 20 
keV), we find 6 max = 0.04 deg . Therefore, the Coulomb trajectory is very closely a straight 
line. 

VI. CONCLUSIONS 

We have developed a general theoretical framework through which the collisional excita- 
tions of collective plasmon resonances in metallic clusters can be quantitatively described. 
The idea of a doorway state representing the resonance, which is coupled to both the open 
decay channels and the chaotic degrees of freedom of the cluster, is used for this purpose. 
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Application of the semiclassical scattering theory in conjunction with the "exit" doorway 
model is made to light Na-clusters using electrons and protons as projectiles. Although 
Pauli exchange effects in the electron scattering are important and they result in non-local 
potentials, we have ignored these effects here for simplicity. The formalism can be applied 
easily to cluster-cluster collisions. Furthermore, it can be generalized to cases involving the 
excitation of multiplasmon states in both spherical and deformed clusters. In the deformed 
cluster case one sees a two-peaked one plasmon resonance 0] and should see a three-peaked 



two-plasmon resonance, a four-peaked three-plasmon resonance etc. ||12j| . These resonances 
can be reached in multistep processes (ground state— >■ one plasmon— > two plasmon—* etc.). 

In our theory only the average cross-section corresponding to the excitation of the door- 
way state was calculated. There should be fluctuation contributions, not considered here. 
Such contributions will be discussed in future work. 
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Figure Captions. 

Figure 1: Collision of a charged projectile (open circle) with metal cluster (solid circle). 
In this illustration, the charge of the cluster has the same sign as that of the projec- 
tile. In this reference frame, the trajectory lies on the xy-plane and the z-axis points 
upwards, away from the paper sheet. The coordinates R and 0, and the distance of 
closest approach, 2a are indicated. 

Figure 2a: Time evolution of the elastic and inelastic probabilities in collisions of 100 
eV electrons with a Na-cluster. The time is indicated in atomic units (a.u.), defined 
as the ratio of the radius of the ground state Bohr orbit in the Hydrogen atom to the 
electron speed in this state. 

Figure 2b: Same as figure 2a for 100 keV protons. 

Figure 3a: The differential cross section as a function of the cluster excitation energy 
in the plasmon region. The projectile is a 100 eV electron. For details see the text. 
The cross section is given in atomic units, which in this case are defined as the radius 
square of the ground state Bohr orbit in Hydrogen. 

Figure 3b: Same as figure 3a for 100 keV protons. 

Figure 4a: Cross sections for the excitation of the plasmon resonance in ionized Sodium 
clusters. The results are given as functions of the collision energy, for a few different 
values of the cluster's net charge. For details see the text. 

Figure 4b: Same as Figure 4a for protons. Here the curves for different cluster charges 
cannot be distinguished. 



18 




Figure 1 




Figure 2a 




time (a.u.) 



Figure 2b 



105 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



100 eV electrons 



10 4 - 



3 io3 



m 10 2 h- 
D 



10 1 - 



Total 



100 




j i i i i i i i i i i i i i i i_u i i 



0.00 



0.05 



0.10 
8 (a.u.) 



0.15 



0.20 



Figure 3a 



103 - 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



100 keV protons 




J I I I I I I I I I I I I I I CjI I L 



0.00 0.05 



0.10 
8 (a.u.) 



0.15 0.20 



Figure 3b 




Figure 4a 




Figure 4b 



